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1 .0 INTRODUCTION 

This report presents C-band and S-band radar error statistics that are recommended 
for use with the groundtracking programs used to process Space Shuttle tracking 
data. The statistics are divided into two parts: bias error statistics, using 

the subscript B, and high-frequency error statistics, using the subscript q. 

Bias errors may be slowly varying to constant. High-frequency random errors 
(noise) are rapidly varying and may or may not be correlated from sample to sample. 

Bias errors are mainly due to hardware defects and to errors in correction for 
atmospheric refraction effects. High-frequency noise is mainly due to hardware 
and due to atmospheric "scintillation". This report identifies three types of 
atmospheric scintillation: horizontal, vertical, and line of sight. This is the 

first time that horizontal and line-of-sight scintillations have been identified. 


2.0 C-BAND ERROR SfATISTICS 


The C-band error statistics apply to the following radars: FPS-16, TPQ- 18 , FPQ-6, 

FPQ-13, FPQ-14, and FPQ-15. Hardware error statistics for the C-band radars for 
the Eastern Test Range were obtained from references 1, 2, and 3. which gave re- 
sults based on 9 months of tracking data. Hardware error statistics for the 
Western Test Range C-band radars were obtained from reference 4, which gave re- 
sults based on 13 months of tracking data. Differences between Eastern Test 
Range and Western Test Range statistics were small. Also, all the various radars 
seemed to have about the same statistics. Note though, that some radars can track 
at longer ranges than others. 


The table below lists the hardware error statistics, subscript H. The symbol 
0q denotes high-frequency error standard deviation. The symbol denotes 

the bias error standard deviation. Subscripts p, A, E denote range, 
azimuth angle, and elevation angle. Maximum and minimum values of a are 
also shown. These are the extreme values observed among the 15 tracking sites 
used to determine the average (RMS) statistics and are generally based on 
several tracking passes. Larger or smaller values may occur for any individual 
pass of tracking data. There also seemed to be little difference between 
the small amount of skin tracking data statistics and the beacon tracking 
statistics. An exception was agp t which seemed to be less. The Ogp for 
beacon tracking was larger, probably due to unknown time delays in the beacon, 
and is the one shown in the following table. 


a qpH = ^*7 m 
OqAH =0.10 mrad 
OqEH = 0.11 mrad 

a BpH = 12 - 5 m 
a bah = 0.08 mrad 
°BEH = mrad 


°qpH,MAX 
°qAH ,MAX 
°qEH,MAX 
a BpH,MAX 
a BAH,MAX 
a BEH,MAX 


4.9 m 
0.14 mrad 
0.15 mrad 
17.4 m 
0.11 mrad 
0.21 mrad 


OqpH.MIN 
°qAH ,MIN 
°qEH,MIN 
°BpH,MIN 
a BAh,MIN 
a BEH,MIN 


s 1.0 m 
= 0.07 mrad 
= 0.07 mrad 
= 7.9 m 
= 0.05 mrad 
= 0.05 mrad 
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The errors listed in the preceding table are the dominar t errors for elevation 
angles above 5 degrees. At low-elevation angles, Oq is increased because of 
atmospheric scintillation noise, rhe magnitude of the scintillation noise is 
dependent on the amount of atmosphere between the tracker and the target. 

The parameter a is used to quantify this amount of atmosphere and is defined 
below. 


Ap = range refraction correction 

Ap * P MEASURED “ P GEOMETRIC 

R 0 = 1 Earth radius = 6 378 165 m 

a = V 1 " exp(-(10 5 Ap/R o ) 4 ) (2) 

A p has an approximate range of values of 
0 < Ap < 144 m 


The parameter a ranges from zero to one. The symbol a = 0 means little 
atmosphere between the tracker and target. The a = 1 means a large amount 
of atmosphere. The a is small, a <.13, for E > 5 degrees. 

To obtain the total a q for the high-frequency noise, the hardware noise is 
root sum squared with the scintillation noise in the following manner. 


Oq = ^Oq H (aO q s) 2 (3) 

Based on personal analysis of low-elevation angle S-band data, the scintillation 
standard deviations are 


<jqps = 0.005 meters 
°qAS = 0.15 mrad 
OqEs = 0.5 mrad 


References 1, 2, and 3 also indicate a value of a q ES 0*5 mrad. Using the 
root mean square (RMS) o q H errors (left column), a table of a q p , °qA and 
OqE versus EmeasURED (tab 16 D has been constructed using a modulus of refrac- 
tion of Nq = 0.000395 and an orbital altitude of 108 n. mi. 


2 
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TABLE I.- 

Oqp, OfqA AND o q E 

FOR H 

= 2-105 METERS 

Em. 

deg 

°qp • 

m 

OqA* 

mrad 

OqE» 

mrad 

0 

2.7 

0.18 

0.51 

.1 

2.7 

.18 

.51 

.2 

2.7 

.18 

.51 

.3 

2.7 

.18 

.51 


2.7 

.18 

.51 

.5 

2.7 

.18 

.51 

.7 

2.7 

.18 

.50 

1 

2.7 

.17 

.47 

2 

2.7 

.13 

.29 

3 

2.7 

.11 

.19 

4 

2.7 

.10 

.15 

5 

2.7 

.10 

.13 

7 

2.7 

.10 

.12 

10 

2.7 

.10 

.11 


Low-elevation angle tracking also causes an increase in the range and elevation 
angle bias error standard deviations. This increase is because of errors in 
correcting for refraction effects. These errors can be quite large. Two 
sources of refraction errors are identified. The first is the error in the esti- 
mate of the modulus of refraction Nq at the observer (station). For Space 
Shuttle tracking, mean monthly values of Nq will be used. A study of refer- 
ence 5 reveals that Nq is changed from month to month by an average amount of 
1.7 percent with a maximum change of 7.5 percent observed at one station. It is 
not meant to imply that using mean monthly averages is bad. In fact, they may 
represent the total atmosphere above the station better than an estimate made 
near the track time. Note that diurnal variations in Nq at inland stations 
may be as much as 6 percent from the mean daily value. A strong weather front 
passage can change Nq by 14 percent. Thus, a 1o error in Nq of 2.5 per- 
cent is suggested for use with Space Shuttle tracking. That is, 


3 
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°N0^^0 s 0.025 

Using the equations of referenoe 6 for low-elevation angles and high altitudes, 
the following approximations can be made. The bias error standard deviation for 
range, due to errors in Nq, is approximately 

°BpN = 0.5(omo/No)Ap (4) 

And for elevation angle 

oben = (o N0 /n 0 )ae (5) 


where Ap and AE are the computed refraction corrections. 

Another major source of bias error due to refraction effects is in the 
algorithms used to compute Ap and AE. These algorithm standard deviations 
are denoted by a Bn A 210(1 obeA- The accuracy of any particular algorithm may 
be checked using tne extensive tables of "exact” corrections in reference 6. 

Even very precise algorithms will contain come error because of the assumptions 
made about the atmosphere, namely that the modulus of refraction N decreases 
exponentially with increasing altitude, and that the Earth is spherical in the 
vicinity of the tracker. The assumption error may be on tlie order of 1 percent 
of the computed correction. The best algorithm for AE in reference 6 had a 
maximum error of 41 percent for E : 0°. At E = 2° the maximum error was 5.85 
percent for low altitudes and 1.6 percent at high altitudes. The best algorithm 
for Ap in reference 6 had an overall accuracy of 1.6 percent. In reference 7, 
the best algorithm for AE had an overall accuracy of 1.7 percent. Reference 
8 gives refraction correction algorithms for altitudes above 54 n. mi. Here, Ap 
had an accuracy of 0.44 percent and AE had an accuracy of 0.55 percent. 

Thus, the total bias error standard deviations are 


OBp 


4 


2 

°BpH 


2 

♦ CFfipN 


2 

♦ OfipA 


(6) 


OBA = OBAH 


(7) 


OBE 




2 2 2 

OBEH + OBEN ♦ OBEA 


( 8 ) 


4 
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°BpH = 12 *5 “ 
o gp = 0.0125 Ap 

°BpA = °* 02 Ap 
a BA = °BAH = 0.08 mrad 
OflEH = 0.12 mrad 
°BEN = 0-025 AE 
°BEA s °.02 AE 


Table II shows the total bias error standard deviations for Nq = 0.000395 and 
H = 108 n. mi. It is seen that range and azimuth bias errors are little 
affected by refraction errors. OgE is strongly affected by refraction errors 
for Em < 5 degrees. 
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TABLE II.- 

°Bp. °BA* °BE 

for H = 2*10 5 

METERS 

E M» 

°Bp» 

°BA* 

°BE» 

deg 

m 

tnrad 

mrad 

0 

12.9 

0.08 

0.63 

.1 

12.9 

.08 

.59 

.2 

12.8 

.08 

.56 

.3 

12.8 

.08 

.53 

.4 

12.7 

.08 

.50 

.5 

12.7 

.08 

.48 

.7 

12.7 

.08 

.44 

1 

12.6 

CO 

o 

• 

.39 

2 

12.6 

.08 

.28 

3 

12.5 

.08 

.23 

4 

12.5 

.08 

.20 

5 

12.5 

.08 

.18 

7 

12.5 

.08 

.15 

10 

12.5 

.08 

.14 

20 

12.5 

.08 

.12 


Little is known about the autocorrelation function of the high-frequency noise 
adding to p , A and E. From reference 9, it appears that the range noise is 
essentially uncorrelated at 10 samples/second. However, the noise adding to azi 
muth and elevation angle is correlated at this sample rate. Exponential correla 
tion is generally assumed with a time constant of t 2 seconds. The 
autocorrelation functions is defined by 


♦ j = EUiCt.j) = EU(ti)e(ti - jAT)) (9) 

EUJ = 0 
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2 2 

Note = E[c|J = Oq is the variance of e*, independent cf the saaple rate. 
Let 


b s e -AT/x (10) 

where AT is the saaple interval and t is called the tiae constant. Siaulated 
exponentially correlated noise is generated in the following manner (ref. 10) 


slnt) = o 

(ID 

E(mni) ~ 0 for i i J 
= 1 for i s J 

(12) 

= be^, ♦ o q Jl - b 2 Hi 

(13) 


It can be shown now that the autocorrelation function of 


e is given by 


= OqbJ = $ 0 EXP(-JAT/t) 


(1m) 


Thus the name exponentially correlated random variable, ECRV. Also note that 
E(e^) =0, it is a zero mean random variable. 

Let us briefly obtain a further insight of the meaning of $j. In the absence of 
any information about ei, the best guess of its value is zero with an error of 


Oq in the guess. 

That is 


£ i = 0 

to q 

(15) 

Now suppose that 
best estimate of 

£i_j is given. 
£i is 

Then it can be shown (ref. 10) that the new 


e i = ($j/to) e i-j - (<t>j/4> 0 ) 2 06) 


For example, consider $j/$o < Then the 1o error is 

r -q ^1 - (4»j/<D 0 ) 2 > O.920 q (17) 


That is, the reduction in the uncertainty is leas than 8 percent. Therefore, 
values of $j/$q less than 9.4 may be considered insignificant in this case. 
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3.0 S-BAND ERROR STATI S TICS 

S-band data consist of four measurements: range p; X angle ax? X angle ay; 

and integrated Doppler data. The subscript D is used for integrated Doppler 
ua^ that is equivalent to a range measurement plus a large bias (constant of 
integration). Integrated Doppler data have units of cycles. The conversion fac 
tor to meters is approximately 15 261 cycles/meter, or 15 cycles/mm. 

All the error statistics were obtained from many passes of satellite tracking at 
various stations. The high frequency error statistics o q will be considered 
first. These errors are due to atmospheric scintillation and to hardware errors 
that are fhe dominant errors for high-elevation angles. OqXH and OqYH were 
obtained by sliding, second-order, polynomial fits to 10-recond data arcs (201 
points) using a total of 2201 points. The midpoint fit value was compared to 
the raw data to obtain the error value. Many passes of 2201 points were used. 

A maximum of 2000 midpoint errors were obtained from these passes, from which 
the autocorrelation function for the X and Y angles were obtained. This 
work was performed under the direction of I. M. Salzburg, head of Orbit Opera- 
tions Section, Operational Orbit Support Branch, Goddard Space Flight Center, 
but has not yet been published. W. Lear reduced the raw autocorrelation func- 
tions to an analytical form, determined the high-frequency range and Doppler sta 
tistics from raw satellite tracking data supplied by Goddard, and determined 
low-elevation angle statistics ( E < 5 degrees). Second-order polynomial fits 
to 51-point samples were U3ed for the angle data to generate error residuals. 
Third-order polynomials were necessary for the range and Doppler data. 

The random noise sigma for the S-band range measurements, obtained by the author 
from Goddard data, is summarized below versus measured elevation angle. 


! 

! Em, degrees 
! 


! 


oqp 


meters 


! 

! 


! ! ! 

! 

! ! 

! 


! 

1.1 ! 

1.5 

! 2.0 ! 2.5 ! 

3.0 ! 3.5 

! 4.0 ! 4.5 

! 5.0 


! 

! 


! ! ! 

! 

! ! 

i 


! 

.33 ! 

.36 

! .43 ! .50 ! 

.43 i .41 

! .41 ! .41 

! .45 

1 

! 

! 


! ! ! 

l 

! ! 

! 



Oqp appeared to be independent of elevation angle; i.e., scintillation effects, 
<jq 0 s* were small. The RMS value of the range noise was 0.43 meters with a maxi- 
mum vaiue 0.69 meters and a minimum value of 0.24 meters. Forty two 51 -point 
samples were used. However, reference 11 indicates that the Space Shuttle S- 
band transponder may be more noisy than indicated above. The hardware theoreti- 
cal value is 10 meters, 3o* Thus, for the Shuttle 


OqpH = 3-3 meters 


For the scintillation error (obtained below from the integrated Doppler data) 


OqpS = 0.005 meters 


t 
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a qp 



(aOqpS ) 2 


= 3-3 meters 


( 18 ) 


where a is given by equation 2. 

Reducing raw integrated Doppler data, the following was obtained by the author. 


l" 


! 

! 


! 


! 


! 


! 


! 


! 


! 


! 

E m , degrees 

! 1.1 

! 

1 .5 

! 

2.0 

! 

2.5 

! 

3.0 

! 

3.5 

! 

4.0 

! 

4.5 

! 5.0 


! 

! 

j 


1 


! 


! 


! 


! 


1 


! 


! 

OqD» mm 

! 5.0 

! 

4.0 

! 

3.2 

! 

3.1 

! 

2.5 

! 

2.4 

1 

2.6 

! 

1.8 

! 1.5 


t 

! 





! 


! 


! 


! 


! 


! 



Thirty eight batches of 51 -point data were used to obtain the values in the 
above table. To convert millimeters to cycles, multiply the above values by 15 
cycles/mm. The above figures suggest a hardware value of OqDH = 1.5 mm; how- 
ever, vibration effects may be larger than this. References 12, 13 and 14 list 
RMS residual values for least-squares orbit determinations. The residuals are 
due to random noise, bias errors, and fit errors and thus represent an upper 
bound for the Doppler noise statistics. The reference values gave OqDH = 8 mm, 
a qDH,MAX = 12.7 mm and 0qDH,MIN = 5 mm. Conversely, five batches in the 
authors analysis showed Oq£>H 1 *9 mm; they were removed from the tabular 
statistics above. Thus, it is suggested 


a qDH = 4 mm = 60 cycles 
OqDS = 5 mm = 75 cycles 


The total standard deviation is thus 


a qD = 



< 6.4 mm = 96 cycles 


(19) 


And, for the S-band station in Hawaii at Kauai, HAW3, it is recommended (ref. 15) 


dqDH = 24 mm = 360 cycles 


For the 10 samples/second data, it was found by the author that the range noise 
had 


r 


4>i/4>o = 0.4 (not strongly correlated) 
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For the integrated Doppler data 

♦ = 0- 2 (essentially uncorrelated) 


Based on Goddard's reduction of the S-band angle data 


°qXH =0.15 rnrad ° q XH,MAX =0.40 mrad 

o q YH = 0-09 mrad O q xH f MAX = °- 16 mrad 


°qXH,MIN = 0-°3 mrad 
°qYH,MIN = O- 04 mrad 


The above figures are based on 29 batches of 2201 raw data points. The angle 
data rate was 20 samples/second. 

Both X and Y, angles are affected by scintillation error. The scintillation 
error is root sum squared with the hardware error in the following manner. 

o qX = yjo qXH ♦ [a(3ax/3E)o q EsJ 2 + i.a(3a x /3A)o q AsJ 2 ( 2 °) 


OqY = \OqYH ♦ (.a(3ax/3E)OqEs) 2 ♦ (a(3ax/3A)aqAs) 2 (21) 

where (ref. 16) 


3a x 


-sin ax/cos aY 


3E 

2 

COS c 

<*X 

sin c ay ♦ sin‘ 

»X 

3ax 





3A 

= cos ay 

sin ax/cos ax 


dniY 

-i 

cos 

a x sin ax 


3E 

2 

cos^ 

a X 

sin 2 ax + sin^ 

<*x 

3cxy 





3 A 

= -sin 

“X 



E = 

arcsin 

(cos ax cos ax) 



30tv 

-8 < <8 for E small 

3E 

= 0 for E = 0 


= 0 for E = 0 


= ±1 for E = 0 
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Because the partials are multiplied by a, which is zero for large E, it may be 
approximated 

da^/d A - dcty/^E = 0 
s 1 

Analysis by the author of low-elevation angle S-band data shows: 


Oq as = 0.15 mrad 
OqES = 0.5 mrad 


Note that at low-elevation angles, Oq/^s represents a horizontal scintillation 
and OqES represents a vertical scintillation, while Oqps and Oqpg repre 
sent scintillation along the line of sight. 

The high-frequency X and Y single random errors are strongly correlated. 
Approximately 30 plots of X and Y angle autocorrelation functions have 
been supplied by Goddard Space Flight Center. These plots have been reduced 
to analytical form as jhown below. For the X angle noise 


4> 


= e - JAT/2 .58 


cos 


/360jAT\° 

\TTT/ 


- 0.3185 sin 


/ 360 jAT\° " 

\TTT/ . 


(22) 


where 4 >q = °qX- For the * angle noise 

♦j/*o ■ e-J^ T/ '- 80 [=03 


- 0.4304 sin 


/36W\° 

\ ‘••87 / J 


(23) 


where 4>o = °qY- Figure 1 shows plots of the X and Y angle autocorrelation 
functions . 


Simulated data having the above autocorrelation functions can be generated by 
(ref. 10) 


e i = - a2 e i-2 ♦ *>1^1 + b 2 n i-1 


(24) 
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where 

E(m) * 0 

E(nini) s 0 for i i J 
=1 for i = J 

$ 3 / 4*0 ~ 

$ 2^0 “ (4>i/$o) 

a 2 = 

$2^0 “ ($i/$o)* 

bl /V^= 1/2 yi - a! ♦ a 2 V 1 - a! - a 2 ♦ 

♦ 1/2 + ai ♦ a2 V 1 + a 1 " a2 - 2<J> 1 /<t>o 

62/ V^o" = 1/2 Vi " a 1 ♦ a 2 V 1 “ a 1 - a 2 ♦ 2<|) i / 4>o 

- 1 /2 Vi ♦ a i + a 2 V^ ■*■ a 1 “ a 2 " 2 $i/$q 

For example, if AT = 0.2 seconds, then the simulated X angle noise is generated 
by 

Ci = 1.7962e i _ 1 - 0.8564e i _ 2 o qX (0. 5156m - 0.5057nj_l) (29) 

Y angle noise is generated by 

£i = 1.7305m-1 - 0.8009£i_2 ♦ o q Y (0.5976m - 0.5846m_l) (30) 

If AT = 0.05 seconds, X angle noise is generated by 

ei = 1.958m-1 - 0.962ei ♦ 0.27 29o qX (m - m- 1 * (3D 

and 4>i/4>o = 0.96, 4>i/4>o = 1 •9584> i _ 1 /4>o - 0 .962<t> i _ 2 /4> 0 . The Y angle noise is 
generated by 


( 26 ) 

(27) 

(28) 
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= 1.94l2e^_i - 0.9^60e^_2 ♦ 0.3240Oqy(ni - Hi-i) (32) 

and 4>i/<t>o s 0.9436, 4>i/$o = 1 •9 1 *12<|) i->-|/4>0 ~ 0.94604>-j_2/$o* 


The recommended high-frequency noise standard deviations for S-band data are 
summarized below. 


°qpH = 3.3 m °qpH ,MAX = 10 m (30) 

°qXH = °- 15 mrad °qXH,MAX = 0<l40 mrad 

OqYH = 0.09 mrad OqYH.MAX = 0.16 mrad 

°qDH = ** mn = 60 cycles 

°qDH,MAX = 12.7 ran = 190 cycles 

°qDH,MIN =0.8 mm = 12 cycles 

OqDH = 24 mm = 360 cycles for HAW 3 


°qpH,MIN 

a qXH,MIN 

°qYH,MIN 


0.24 m 
0.03 mrad 
0.04 mrad 


Oqps = 0.005 m 
Oq as = 0.15 mrad 
°qES = °-5 mrad 
a q DS = 5 mm = 75 cycles 

The bias error statistics due to hardware were obtained from references 13. 14, 
15 , 16 and 17 . 


OfipH = 28 m 
OBXH =0.5 mrad 
a BYH = 0.25 mrad 
OfiDH = 7.5 Hz 

= 0.5 mm/ sec 


°BPH,MAX = 61 m 
a BXH,MAX = 1 -®5 mrad 3 
a BYH,MAX =1-66 mrad 3 

°BDH f MIN = 10 Hz 

= 0:67 mm/ sec 


a BpH,MIN = 18 m 
a BXH,MIN =0.10 mrad 3 
°BYH,MIN = 0.19 mrad 3 
a BDH,MIN = 5 Hz 

=0.33 mm/sec 


Reference 16 states that the Doppler (rate) bias error should be modeled as an 
exponentially correlated random variable adding to the Doppler constant of inte- 
gration. The time constant is less than 400 seconds. 


Refraction correction errors increase the hardware bia3 standard deviations at 
low-elevation angles in the following way. 


y Bp 


-4 


2 

°BpH 


2 

+ °BpN 


2 

°BpA 


(33) 


OfipN is given by equation 4 and is the error due to the error in the modulus 
of refraction Nq at the observer (station), is due to algorithm 

errors in computing the refraction correction Ap. Both and 

should be smail compared to the recommended value for ajjpH = 28 meters. 


3 For an uncalibrated station. 
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The total X and Y angle bias standard deviations are 


a BX 


°BY 



( 3*0 


(35) 


°BEN is the error due to errors in Nq and is given by equation 5. OggA 
is due to the algorithm error in computing AE and is a function of the 
algorithm used. 
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Figure 1.- Normalized auto correlation functions. 
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